New method to evaluate divergent series via the Wigner function 
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It is shown how a physical function, namely the Wigner function, that in principle may be mea- 
sured, can be used to evaluate divergent series. 
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We give a physical way to evaluate divergent series. 
Abel quoted on them: 

" Divergent series are on the whole devil's work, and it 
is a shame that one dares to found any proof on them. 
One can get out of them what one wants if one uses them, 
and it is they which have made so much unhappiness and 
so many paradoxes" 

We can evaluate the alternating series S = 1 — 1 + 1 — 
1 + 1 + . . . by considering the alternating geometric series 



1 — x + x 



— ... provided \x\ < 1, we know that 



it converges to 1/(1 + x). By allowing x to get close to 
1, the series will tend to a value close to 0.5. Therefore, 
we can say that the (divergent) sum 5* evaluates 1/2 in 
the Abel sense [l|. 

Another way of defining convergence is by using 
Cesaro's sums. The sum S diverges not because the par- 
tial sums grow uncontrollably, but rather because the 
partial sums oscillate. If we could find a way of averag- 
ing the sums in order to smooth them out, maybe this 
series will converge. Using Abel and Cesaro's limits it 
has been shown also that 1 — 2 + 3 — 4 + 5 + ... has dis- 
crete sum —1/4 0, values that agree with evaluations of 
the Riemann zeta function (2j. 

In this contribution we propose a new method to evalu- 
ate divergent series by using a physical function, namely 
the Wigner function [3J, which may be written in the 
form, [4[ 



W(a) = -Tr{(-lfD*(a)pD(a)} 



(1) 



where p is the system's density matrix, the number op- 
erator n = a^a with a the annihilation operator and 
D(a) — exp(aa^ — a* a) the Glauber displacement op- 
erator, where a = (q + ip)/y/2. Another common form 
for the Wigner function is 



W A (a) = ±- [due m P(q 
2ir J 



■u/2\A\q-u/2) (2) 



where we define it for an arbitrary operator A. If wc 
consider a function of the position operator, we find 



W m (a) 
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due lup f(q - u/2)8(u) 



M 

2tt ' 



(3) 



It is worth noting that the Wigner function is in gen- 
eral a physical function and that has been measured in 
experiments, in particular for the first excited state of 
the vibrational motion of an ion Q . 



For an arbitrary operator, for instance, position to an 
arbitrary power q k the Wigner function may also be given 
in a series representation 



f- = i£(-l)'>|£W£(a)|n) 

71 = 

= i£(-l) n <n|(9 + ?)*!»>• ( 4 ) 

where \n) is a number state. Note that equation (|4|) 
has already the alternating form of the series considered 
above. 

By doing k = 0, we obtain 
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doing k = 2 gives @ 
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J2(-lT(q 2 + (n\q 2 \n)), 
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(5) 



(6) 



which by using that (n\q 2 \n) = n+1/2 allows us to eval- 
uate the sum 



E(- 1 ) B » = -4' 

71 = 



(7) 



where we have used ([5]) to evaluate ([Tj) . Eq. (J6J) then 
shows that there is a recursion relation between the 
higher order sums in terms of the lower order ones. We 
can find a general expression that will have this recur- 
sion in it. In order to do this we write q in terms of 
annihilation and creation operators, 



<l 



a + a 1 

"7T' 



= y/n\n-l) 



'» = sjn + l|n + 1), 



(8) 



and insert it in (J4J) , this gives 



^ = £(-D n £ 



71=0 



s=0 



k \ q 



k—s 



2 s/2 



{n\{a + a?) s \n). (9) 



by equating the coefficients of powers of q at right and 
left of the equal sign, we finally obtain (s > 0) 



= E(- 1 )"H(« + « t ) 2s |n). 



(10) 



n=0 
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This equation is the main result of the manuscript, as it 
will allow the evaluation of the divergent series considered 
earlier. 

Because of the average with the numbers states, we can 
remove all terms that do not contain an equal number of 



tribute to the sum of diagonal matrix elements. There- 
fore by considering those elements we may neglect all 
other terms of the sum. We have Q 



(a 1 + a y 



( m ™ ) : (a t o) W2 :w> m 







odd 



(11) 

where : (a^a) s :w denotes the Weyl (symmetric) ordering 
of the operator h s . It may be transformed into normal 
ordering using [|[ 



2=0 



(a'aj : w = Z^^i [ i I a a • ( 12 ) 



Inserting (fT2|) into (| 10[) and making use of the expression 



«*V» = ^-^|n>, (13) 



s and a's, as they are the only terms that will con- we obtain 



EE 

n=0 1=0 



(-iy 



2 l \l \s-l 



= 0. 



(14) 



For s — 1 we obtain ([5]) and for s = 2 we obtain 0, 
while for s > 2 we obtain the sums 



(15) 



n=0 



in terms of the lower sums, i.e. as a recursion formula. 

In conclusion, we have shown a physical form, to eval- 
uate some divergent sums, by means of a function that in 
principle may be measured, namely, the Wigner function. 
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